Loops were firstly called using the ABC Model. Since this method integrates ATAC-seq, H3K27ac and gene expression data when calling enhancer-promoter (E-P) pairs, we can be more confident we are truly calling loops. After defining E-P loops, differential loop calling was performed using a modified DESeq2 approach via Hictoolsr.

Activity (A) is calculated as the geometric mean of the read counts of ATAC-seq and H3K27ac chromatin immunoprecipitation sequencing (ChIP–seq) at element E, and Contact (C) as the KR-normalized Hi-C contact frequency between E and the promoter of gene G at 5-kb resolution
Download Enhancer_promoter_loops.csvDifferential loop analysis was performed at 1000kb resolution (ie. regions were binned at 10kb windows).

